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Abstract. We report on experimental studies of divergence of proton beams from 
nanometer thick diamond-like carbon (DLC) foils irradiated by an intense laser with 
high contrast. Proton beams with extremely small divergence (half angle) of 2° are 
observed in addition with a remarkably well-collimated feature over the whole energy 
range, showing one order of magnitude reduction of the divergence angle in comparison 
to the results from /im thick targets. We demonstrate that this reduction arises from a 
steep longitudinal electron density gradient and an exponentially decaying transverse 
profile at the rear side of the ultrathin foils. Agreements are found both in an analytical 
model and in particle-in-cell simulations. Those novel features make nm foils an 
attractive alternative for high flux experiments relevant for fundamental research in 
nuclear and warm dense matter physics. 
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1. Introduction 

The emission of highly energetic ions from solid targets irradiated by intense laser pulses 
has attracted great attention over the past decades [1] . The short time scale on which the 
acceleration occurs along with the small source size enable extremely high ion densities 
in the MeV bunches which could be superior for specific applications [2J. However, 
such high density is only maintained close to the source, it drops quickly due to the 
angular spread of few tens of degrees [3j Hj. Such large angles lead to large losses 
using magnetic quadrupoles (5], complicate the beam transport and therefore trigger 
investigations on sophisticated transportation schemes such as pulsed solenoid [6] and 
laser driven micro lenses [7]. Meanwhile, shaped lens target [8], droplets [9] and curved 
target [10] have been used to manipulate the ions angular distribution. Those approaches 
were mostly based on target normal sheath acceleration (TNSA) [TT] with /im thick 
targets. Acceleration fields are built at the target rear by the hot electrons generated at 
the front side of the targets. The divergence of the ions strongly depends on the electron 
density and phase space distribution of the electrons behind the target, which is initially 
related to the laser profile and then disturbed during the transportation through the 
targets [12] . 

Recently, ultrathin foils with thickness down to nm scale have been investigated 
experimentally p3J [HI EE] > enabled by the improvement on the laser temporal contrast. 
A divergent ion beam with opening half angle of roughly 10° from 50 nm ultrathin foils 
can be inferred from [Hj. Proton beams with divergence of 5-6° have been observed 
from 800 nm CH targets [15]. These results imply that thinner foils can generate much 
more collimated ion emission as compared to fim thick targets. 

In this letter, we present the first detail study of the divergence of proton beams 
accelerated from 5-20 nm thick Diamond-Like-Carbon (DLC) foils [16]. Divergences 
as low as 2° were observed for different foil thickness and irradiation conditions. The 
proton beams show a pronounced collimation over the whole energy range. We attribute 
the small divergence to a steep longitudinal electron density gradient at the target 
rear in comparison with /im thick target, and the collimation feature is related to an 
exponentially decaying electron density in transverse dimension. Our interpretation is 
supported by both an analytical model and the two-dimensional particle-in-cell (PIC) 
simulations, revealing new physics that occur in the interaction of high-intensity laser 
with nanometer thin targets. 

2. Experiment 

The experiments were performed with the ATLAS Ti: sapphire laser system at Max- 
Planck-Institute for Quantum Optics. This system delivers pulses with a duration of 
30 fs FWHM centered at 795 nm wavelength. The initial laser contrast is 3 x 10~ 6 at 
2 ps before the peak of main pulse. A re-collimating double plasma mirror system was 
introduced to further enhance the value to 10~ 9 . 400 mJ laser energy was delivered on 
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target. A 90° f/2 off-axis parabolic mirror focuses the pulses to a measured FWHM 
diameter spot size of 3 /im, yielding peak intensity of 8 x 10 19 W/cm 2 . DLC foils of 
thickness 5, 10 and 20 nm have been irradiated under normal incidence for varying spot 
size, the actual spot size on target in the range of 3-19.2 /xm has been adjusted by 
moving the target along the laser axis. 
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Figure 1. (color online). Experimental setup. The divergence of protons characterized 
with the magnetic spectrometer and Image plate (IP). The trajectory of protons 
through the magnetic field as well as the magnetic field structure is shown in the 
setup picture. The resulting isoenergy contours of the magnetic spectrometer are 
superimposed with a raw image of proton energy- angular distribution as it appears at 
the IP. The example is obtained with a 10 nm target displaced by 100 fim from the 
laser focal plane. 

Fig. [1] shows the magnetic spectrometer with a gap of 14 cm employed for the 
proton beam measurements. A long vertical entrance slit of 300 fim width is placed in 
front of the magnetic field. This configuration enables angularly-resolved high-accuracy 
energy distribution measurement. Fujifilm BAS-TR image plates (IPs) were positioned 
at a distance of 30 cm behind the magnets to capture ion phase space over an angular 
range of 8°. The IPs have been absolutely calibrated at MLL Tandem accelerator [17J. 
A Layer of 45 /im Al foil was added in front of the IPs to block heavy ions and to 
protect IPs from direct and scattered laser light. Protons with energies beyond 2 MeV 
are recorded. A typical raw image from a 10 nm DLC foil is shown in Fig. [TJ The image 
of the entrance slit, i.e., the zero line, and the low energy cutoff line from proton signal 
allows to extrapolate the average magnetic field for different angles. These are used 
to transfer the spatial information from raw image to the energy angular distribution 
of protons. The trajectory of proton beams through the dipole magnets as well as the 
magnetic field was carefully modeled and calibrated, as shown in Fig. [TJ The resulted 
isoenergy contours of the spectrometer shows how the proton propagate through the 
spectrometer and form a angular distribution on the detector, the isoenergy curve shows 
fair agreement of the experimental raw image for 2 MeV low energy cutoff by the Al 
foil. 

The smallest divergence was observed with a 10 nm target displaced by 100 /xm 
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Figure 2. (color online), (a) Beam divergence half angle as a function of proton 
energy for the most collimated beam from DLC foils, along with other data published 
in the literature from /im thick targets (the light blue area, presented by blue, green, 
black and cyan curves) [21 (HI H] and from nm scale thin targets (brown dot and magenta 
curve) [21 [15]. (b) Experimentally processed result of the data presented in Fig. [2ja) 
after normalization, where the color scale denotes the normalized fluence for given 
energy. 

from the laser focal plane. Fig. Mjo) shows the energy-angular distribution of the 
example from Fig. [1] after normalization in order to highlight the collimation over the 
complete energy range. The divergence is almost constant over the detected energy 
range, resembling the afore mentioned collimation. The angular distribution is fitted by 
a Gaussian function for each energy value. We define our divergence by the half value 
of Full-width half-maximum (FWHM) of the fitting profile and plot these values as a 
function of proton energy normalized to the peak energy, as shown in Fig. [2]^a) by the 
red curve. This enables a comparison to established results [3j EJ E] that also presented 
in Fig. |2]^a). Compared to jum thick targets, the half angle from 10 nm foil is reduced 
by a factor of 10. Moreover, the typical increasing of divergence with decreasing energy 
is not observed. The results from experiment with 50 nm [TJ] and 800 nm foils are 
included as well, indicating an overall reduction of divergence with decreasing target 
thickness. 
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Figure 3. (color online). (a)Beam divergence (half angle) for varying thickness of 
DLC foils and target positions. The lower axis shows the target position, where +/- 
means the foils were placed before/after laser focal plane, as indicated by the small 
pictures. While the upper axis denotes the focal spot FWHM diameter given by 
the spot size of a perfect Gaussian beam in the propagation axis and quantified by 
the measurements from real laser intensity distribution. Here the vertical error bar 
indicates the standard deviation for each shot while the horizontal error bar shows 
the positioning accuracy of about 10 /jm, smaller than the Rayleigh length of 25 (im 
marked by the light blue area. The black curve is an empirical fitting curve showing 
the divergence scale with laser FWHM diameter (Dl)^ 1 / 2 . The inset presents the 
corresponding total number of protons for each shot, extracted from the measured 
proton spectrum above 2 MeV and corresponding divergence from the parameter scan. 
The Monitored laser intensity distribution is shown at three target positions -50 /im 
(b), fim (c) and 50 ^m (d). 



Further on, the thickness of the targets and their positions with respect to the focal 
plane of the laser were varied. As shown above, the divergence showed no noticeable 
dependence on energy. Therefore, we plotted the average value of the half angle as 
a function of the target position in Fig. |3]^a). The divergence is maximized with a 
value of 4.6° in the focus plane for the thickest foil of 20 nm. For thinner foils (5 and 
10 nm), these values are reduced to 3.3°. The divergence of the protons decreases with 
increasing focal spot size on the target when moving the target out of focal plane to both 
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sides beyond the Rayleigh length. The smallest divergence was obtained in the target 
position of +100 fim, which is our example of Fig. 12(a), (b). Additionally, the total 
number of protons above 2 MeV even increases when enlarging the target positions, i.e. 
when the divergence reduces, as shown in the inset of Fig. [3(a). Note that the spectral 
distributions remained exponential throughout the whole parameter scan, in contrast 
to the case of the previous work with circular polarization [131 EE] or under much higher 
laser intensity with linear polarization from experiments [19] . 

3. Discussion and outlook 

To explain our results, we derive the emission angle of ions 9 = arctan(/ °° E y dt/ J °° E z dt) 
based on a quasi- stationary electrostatic model, where y and z are the transversal 
and the longitudinal dimension, respectively. The electric field strengths are deter- 
mined by E = — V$. With a Boltzmann distribution of the hot electrons, the elec- 
tric field is then given by E y oc — ^p^, E z oc -^^^§7, where n e is the elec- 
tron density. Defining the local field direction as a = E y /E z , the emission angle 
9 = arctan(/ °° aE z dt/J °° E z dt) « arctan<a>= arctan<^/^j>, where the angle 
bracket denotes the average along the ion trajectory. We assume the electron density 
n e — no ■ £(y) ■ exp[— f-] with a transverse profile £(y) and a longitudinal exponential 
distribution, as widely accepted [20] . Here n is a constant denoting the electron density, 
and l is the longitudinal density scale length. One derives 

am 

9 = arctan < — ■ l >, (1) 

The emission angle 9 represents the characteristic value of the measured divergence. 

Eq. (Op) shows how the divergence relies on the electron density distribution. It not 
only predicts the influence of transverse electron density on the divergence of ions, 
consistent with [12], but highlights the importance of longitudinal electron density 
distribution as well. This fact has not been well clarified so far. A steep longitudinal 
electron density gradient, i.e., small Iq, will reduce the divergence of ions. 

The scale length Iq can be estimated as the local Debye length Ad = 
(eo/csT , e /'^oe 2 ) 1 ^ 2 [2I], where ksT e is the electron mean energy. For a given fcsT e , l 
is inversely proportional to the square root of the electron density n . In the case 
of fim thick targets, due to the large angular spread of few 10s of degrees for hot 
electrons inside the target [22], no drops significantly at the target rear, resulting in 
a typical Iq of few /im [20]. When the target thickness is reduced from tens of fim to 
submicrometer, the influence from electron propagation through the target is suppressed 
and the angular spread of hot electrons is substantially reduced, altogether resulting in 
a higher no- Moreover, with the reduction of thickness, the recirculation of the hot 
electrons is enhanced [23], which can further increase n$. Those arguments imply that 
a small Iq exists at the rear side of thinner targets, which in turn resulting in a small 
divergence of ions. When the thickness of the target is down to nm scale, an additional 
factor arises; The pondermotive force pushes a large portion of electrons away from the 
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Figure 4. (color online), (a) Transverse electron density distribution (y axis) at 
three different z positions (as denoted by different color lines in Fig. 4(b)) and (b) 
Longitudinal electron density distribution (z axis) at y = at t = 30T (when the 
laser pulse has just left the target), where T corresponds to one laser cycle. Those 
distributions are obtained with 40 nm target irradiated at the best focus position. 
The front side of the target was initially located at z = 10 fim. The dots denote the 
simulated results while the curves are the exponential functions, (c) Simulation result 
for proton angular distribution at t = 100T after normalization as a comparison to 
the experiment result in Fig. Hfb). Here the red rectangular marks the experimental 
observation window. 

ions, which further increase the electron density as well as the density gradient, in turn 
reducing l . Ultimately in light sail regime, a balance between ponderomotive force and 
charge separation field is formed, and a dense electron layer pushed by the laser pulse 
drives the ion acceleration, the acceleration field is confined between two oppositely 
charged plates within 10s of nm [2H 125] . 

In previous experiments, the transverse electron density was found to be well 

2 

approximated by a bell shape profile £(y) oc exp(— 4 In 2-^) for fim thick target [121 [27] . 
where D corresponds to the FWHM diameter of the transverse Gaussian profile. Note 
that typically D = Dl + 2d tan {3 is assumed [26] , which is much bigger than the FWHM 
diameter of laser spot Dl due to the fim target thickness d and large half angle of the 
electron angular spread (3. In this case Eq. ([T]) reads 8 cx arctan<81n2-^- • l >. It is 
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monotonically increasing with local transverse position y. Since the high acceleration 
field appears in the centre, the high-energy ions originate from small y. This implies 
a reduction of divergence with increasing energy, consistent with [31 El HI [28] but 
contradicts our results. Regarding Eq. ([1]), the collimation manner observed in 
our experiments can be explained by an exponentially decaying transverse density 
distribution £(?/) cx exp[— where l y is the transverse density scale length. With 
such a given transverse profile, a constant value of 6 cx arctan(Z // y ), independent on the 
energy and the trajectory of ions, is obtained, which is consistent with our observations. 
For ultrathin targets, much more low-energetic electrons can penetrate through the nm 
foil and contribute to the acceleration field. Such low-energy electrons are trapped close 
to target and reach equilibrium in a short time and form a Boltzmann distribution along 
the transverse dimension as well, leading to a constant divergence as observed in the 
experiment. 

To verify our arguments and analytical model, 2D PIC simulation were performed 
with the KLAP code [25]. Solid density (no = 350n c , where n c is the critical density) 
plasma slab of 40 nm thickness was considered. The initial temperature of electrons is 
1 keV. The simulation box is 60 A in laser direction (x) and 20 A in transverse direction 
(y) in 2D with a resolution of 200 cells/A and 40 cells/A, respectively. Each cell is filled 
with 400 quasiparticles. A linearly polarized laser pulse with a Gaussian envelope in 
both the spatial and temporal distribution with a FWHM diameter Dl of 3 /im and a 
FWHM duration of 33 fs, is used to replicate the experiment conditions. 

Fig. H(a) shows the transverse electron density profile at three different z positions 
behind a 20 nm foil. Indeed, the profiles are found to be well represented by an 
exponentially decaying profile with a scale length l y , similar to the radius of laser spot 
D L /2 of approximately 1.5 /im. [see Fig. Eta)]. The longitudinal density distribution 
from the simulation in Fig. H^b) is represented by an exponential function as well. The 
density scale length Iq =1/12 /im, is much smaller than the typical scale length Iq for /j,m 
targets. Such an electron density distribution finally leads to a well collimated proton 
beam [see Fig. E(c)], in a good agreement with the experimental observation in Fig. 
12(b). The divergence does not depend on energy, with a constant value of about 2° over 
the whole energy range. By using l y and Iq from simulation, Eq. (CQ) gives a half angle 
of 3.2°. 

Furthermore, Eq. ([I]) predicts the reduction of 9 with increasing the transverse 
scale length of electron distribution l y , which in the case of ultrathin foils, is related to 
the laser FWHM diameter Dl, as confirmed by our target position scan[cf. Fig. [3(a)] . 
Although, a precise determination of the values of l y and Iq accordingly to different 
-D^ is beyond the capacity of our simple model and our laser intensity distribution 
outside the focal plane is not perfect Gaussian, we found that the divergence half angle 
roughly scales with (Dl)^ 1 ^ 2 experimentally, i.e., Iq/I v oc (Dl)^ 1 ^ 2 - Note that the laser 
intensity distribution was monitored for different target positions during the experiment 
campaign. The laser intensity distribution is close to an uniform Gaussian profile for 
most of the target positions, as shown in Fig. [3(a), (b) and (d). We therefore ignore the 
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influence of possible nonuniformity of the laser intensity distribution on the divergence. 
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Figure 5. (color online), (a) Electron spectrum at t = 30T and (b) Corresponding 
longitudinal electron density distribution (z axis) at y — at t — 30T (when the laser 
pulse has just left the target) with two different laser intensities ao = 20 (red curve) 
and ao = 5 (black curve), where T corresponds to one laser cycle. The front side of 
the target was initially located at z — 10 /im. (c)Simulation result for proton angular 
distribution at t = 100T with ao = 20. 

Our experiment were carried out with a small laser system. To investigate the effect 
of larger laser intensity, further simulations were carried out with identical parameters 
except for a large ao = 20 as compared to our experimental conditions with a = 5, 
where the laser intensity Iq is 16 times bigger as ao is 4 times bigger. The simulation 
results are shown in Fig. [51 the electron temperature is increased by more than 3 (from 
0.73 MeV to 2.4 MeV) with higher laser intensity (Iq is 16 times bigger). Our model 
predicts that the divergence depends on l , while l increases roughly with the square 
root of electron temperature, i.e., from 1/12 /im to 1/7 /im in our example. Therefore we 
expect a two times larger divergence angle, consistent with our simple model. Still, the 
ion beams exhibit an almost constant divergence angle under higher laser intensity. Note 
that the result of proton angular distribution from simulation in Fig. |5] was normalized 
to the maximum proton energy, where the maximum proton energy is increased from 
6 MeV to roughly 30 MeV by a factor of 5 with higher laser intensity. Our simulation 
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indicates that the divergence depends only weakly on the laser intensity with other 
parameters unchanged, i.e., 9 oc (Jo) 1 ^ 4 - 

4. Conclusion 

In conclusion, we investigated the divergence of proton beams generated from ultrathin 
DLC foils in detail. We demonstrated well collimated proton beams with divergence half 
angle as low as 2°. This constitutes the smallest value reported so far and one order of 
magnitude lower than achieved with jum targets. As a consequence, 100 times increase 
in proton fluence is observed [29]. Moreover, the proton beams are well-collimated 
over the complete energy range and can be further optimized by adjusting the focal 
spot size. Those ultrasmall divergence and well collimation is considered as an intrinsic 
characteristic of nm ultrathin foils. We expect these observations to be of particularly 
interest for applications. For example, the divergence angle is sufficiently small to 
enable lossless beam transport utilizing conventional ion optics such as quadrupoles with 
typically a few msr acceptance angle [5]. More directly, investigations that rely on high 
proton flux, can benefit from the ultrasmall divergence, possibly enhanced by focusing 
laser accelerated ions with shaped nm targets |30j , a development that is pursued in our 
laboratory. 
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